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ABSTRACT 


The  Defence  Research  Establishment  Valcartier  currently  supports  the  develop¬ 
ment  of  a  generic  six-degree-of-freedom  simulation  model  of  a  tactical  air-to-air  missile.  For 
efficiency,  the  development  of  this  missile  simulation  takes  advantage  of  already  existing 
model  components  when  available,  one  example  being  the  modeling  of  the  missile  dynam¬ 
ics.  The  documented  six-degree-of-freedom  missile  dynamics  model  currently  available  uses 
a  non-standard  system  of  axes,  where  the  x  axis  is  forward,  the  z  axis  is  lateral  and  the  y 
axis  is  downward,  to  describe  the  missile  body,  seeker  and  missile-target  line-of-sight  vector. 
However,  the  most  common  axes  system  in  current  simulation  models  and  in  the  DREV 
library  of  weapon  system  model  components  is  the  NED  system  (North-East-Down),  which 
uses  a  forward  x  axis,  a  downward  z  axis  and  a  lateral  y  axis.  In  order  to  maintain  com¬ 
patibility  with  the  existing  library  of  components,  this  missile  dynamics  model  has  been 
modified  to  use  the  NED  standard  system.  The  mathematical  relations  of  the  resulting 
model  are  presented. 


RESUME 


Le  Centre  de  recherches  porn-  la  defense  de  Valcartier  soutient  pr4sentement  le 
developpement  d’un  module  g^nerique  &  six  degres  de  liberte  d’un  missile  tactique  air- 
air.  Dans  un  but  d’^fficacite,  le  developpement  de  ce  modele  de  missile  prend  avantage  de 
composantes  d6ja  existantes,  comme  par  exemple,  dans  le  cas  de  la  dynamique  du  missile. 
Le  module  de  reference  &  six  degres  de  liberte  de  la  dynamique  du  missile  tactique  utilise 
presentement  un  systeme  d’axes  different  du  standard,  soit  1’axe  x  vers  Pavant,  l’axe  y  vers 
le  bas  et  l’axe  z  en  lateral,  pour  representer  le  corps  du  missile,  l’autodirecteur  ainsi  que 
la  ligne  missile-but.  Le  systeme  d’axes  le  plus  couramment  utilise  dans  le  domaine  de  la 
simulation  de  missiles  tactiques  ainsi  que  dans  la  librairie  de  composantes  de  module  de 
simulation  du  CRDV  est  le  systeme  d’axes  NED,  qui  utilise  1’axe  x  vers  1’avant,  l’axe  z 
vers  le  bas  et  l’axe  y  en  lateral.  Pour  assurer  une  complete  compatibilite  entre  le  module 
d6velopp6  et  la  librairie  de  composantes  existante,  celui-ci  a  ete  modifie  pom:  utiliser  le 
systeme  d’axes  standard  NED.  Le  developpement  des  relations  matbematiques  du  modele 
modifie  est  presente. 
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EXECUTIVE  SUMMARY 


A  three-degree-of-freedom  computer  simulation  model  is  currently  used  within 
DND  to  assess  the  performance  of  tactical  air-to-air  missiles.  This  simulation  capability 
is  not  of  a  sufficiently  high  fidelity  to  support  CF-18  missiles  operational  capability,  and 
a  requirement  for  a  high-fidelity  simulation  model  has  been  defined.  In  response  to  this 
requirement,  high  fidelity  and  validated  six-degree-of-freedom  missile  simulation  models  are 
being  developed  to  accurately  model  the  aerodynamics,  guidance  and  control  characteristics 
of  current  tactical  air-to-air  missiles  for  system  performance  assessment. 

For  efficiency,  the  development  of  this  missile  simulation  takes  advantage  of  already 
available  model  components,  one  example  being  the  modeling  of  the  missile  dynamics. 
However,  the  documented  six-degree-of-freedom  model  of  the  missile  dynamics  currently 
available  is  not  based  on  the  standard  NED  (North-East-Down)  axes  system  used  in  most 
current  simulation  models  and,  in  particular,  in  the  DREV  library  of  weapon  system  model 
components.  In  order  to  maintain  compatibility  with  the  existing  library  of  components, 
this  missile  dynamics  model  has  been  modified  to  use  the  NED  system.  The  definition  of 
the  missile  dynamics  model  in  the  NED  system  represents  a  preliminary  milestone  towards 
the  acquisition  of  the  required  simulation  capability. 
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NOMENCLATURE 


axbi  aybi  azb 

- 

missile  body  acceleration  components  (BF) 

BF 

- 

missile  body-fixed  coordinate  system 

dla 

- 

distance  from  the  center  of  gravity  to  the  body  accelerometer 

dmiss 

- 

miss  distance 

DND 

- 

Department  of  National  Defence,  Canada 

DOF 

- 

degree-of-freedom 

DREV 

- 

Defence  Research  Establishment  Valcartier 

EF 

- 

earth-fixed  coordinate  system 

Fx,  Fy,  Fz 

- 

sum  of  external  forces  acting  on  the  missile  body 

g 

- 

gravitational  constant 

hj' 

- 

target  altitude 

H{s) 

- 

missile  seeker  tracking  loop  transfer  function 

1 xxi  lyyi  Izz 

- 

moments  of  inertia  of  the  missile  body 

Kp 

- 

tracking  loop  gain 

LF 

- 

LOS-fixed  coordinate  system 

LOS 

- 

missile-to-target  line-of-sight 

Mm 

- 

missile  mass 

Adxi  Ady,  z 

- 

sum  of  aerodynamic  moments  acting  on  the  missile  body 

NED 

- 

NED  system  of  axes 

R 

- 

Range  vector 

F8a[]iRBe[ 

- 

Euler  rotation  matrix 

SF 

- 

missile  seeker-fixed  coordinate  system 

tut2 

- 

transformation  matrices 

TF 

- 

target-fixed  coordinate  system 

U,  V,  w 

- 

components  of  the  missile  velocity,  Vm  (BF) 

Vm 

- 

missile  velocity 

Vt 

- 

target  velocity 

VF 

- 

missile  velocity-fixed  coordinate  system 

Xi,Vi,Zi 

- 

cartesian  axes  of  coordinate  system  i 

®i ?  Vii  Zi 

- 

unit  vectors  along  the  cartesian  axes  of  coordinate  system  i 

&l\±b 

- 

contribution  of  Xb  to  ii 

xp 

- 

unit  vector  perpendicular  to  the  plane  formed  by  xs  and  xi 

Subscripts 


A  -  azimuth 

b  -  missile  body-fixed  coordinate  system 

E  -  elevation 

il  -  intermediate  I  coordinate  system 


i2 

1 

M 

MT 

MT-1 

r,  R 

RL 

s 

T 

t 

UD 

v 


Greek  Symbols 


C 

&A' )  aE 

CXRL,CtUD 

—* 

e 

£RL,£UD 

$ 

<pA,<t>E,<i>R 

$1 

4>r,  <!>UD,  <j>RL 

<f>RL,  <f>UD 
'YA,'YE,7R 

rA,rE,rR 

A RL  J  A UD 

0A,0E,OR 

Qti9vd,0rl 

<7A,0E,<7R 

<?r,<?UD,GRL 

tib 

UJl 

&l/s 

d>, 

W./6 

U>xei  ^yei  ^ze 
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-  intermediate  II  coordinate  system 

-  LOS-fixed  coordinate  system 

-  missile 

-  missile-target  (EF) 

-  missile- target  (LF) 

-  roll 

-  lateral  (right-left) 

-  missile  seeker-fixed  coordinate  system 

-  target 

-  target  body-fixed  coordinate  system 

-  vertical  (up-down) 

-  missile  velocity-fixed  coordinate  system 


-  projection  of  a#  on  the  earth-fixed  xy  plane 

-  Euler  angles  of  the  missile  body  intermediate  II  coordinate  system  (VF) 

-  Euler  angles  of  the  missile  body  (VF) 

total  angular  error  vector  from  the  seeker  boresight  to  the  LOS  vector 
Euler  angles  of  the  LOS  vector  (SF) 

seeker  tracking  rate  vector 

-  Euler  angles  of  the  missile  seeker  (EF) 

-  limited  seeker  tracking  rate  vector 

-  missile  seeker  angles  (SF) 

-  components  of  < p  (SF) 

-  Euler  angles  of  the  missile  velocity  vector  (EF) 

-  Euler  angles  of  the  target  velocity  vector  (EF) 

-  Euler  angles  of  the  seeker  (BF) 

-  Euler  angles  of  the  missile  body  (EF) 

-  missile  body  angles  (BF) 

-  Euler  angles  of  the  LOS  (EF) 

-  LOS  angles  (LF) 

-  missile  body  angular  rate  vector 

-  LOS  angular  rate  vector 

-  LOS  angular  rate  vector  relative  to  the  seeker 

-  missile  seeker  angular  rate  vector 

-  missile  seeker  angular  rate  vector  relative  to  the  missile  body 
components  of  the  LOS  angular  rate  vector  (EF) 
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1.0  INTRODUCTION 

To  effectively  contribute  to  improving  the  CF-18  missiles  operational  capability, 
the  existing  simulation  model  and  associated  modeling  environment  in  use  within  DND  to 
evaluate  air-to-air  missile  system  performance  require  improvement  with  regard  to  fidelity, 
accuracy  and  reliability.  The  field  of  operational  research  and  analysis  would  also  benefit 
from  improved  models  in  the  area  of  performance  trade-off  studies  (options  analysis)  to 
support  the  acquisition  of  new  missile  systems.  In  response  to  this  deficiency,  a  high  fidelity 
and  validated  six-degree-of-freedom  missile  simulation  model  is  currently  under  develop¬ 
ment  to  accurately  model  the  aerodynamics,  guidance  and  control  characteristics  of  current 
tactical  air-to-air  missile  systems. 

For  efficiency,  the  development  of  this  missile  simulation  takes  advantage  of  al¬ 
ready  existing  model  components  when  available,  one  example  being  the  modeling  of  the 
missile  dynamics.  The  available  documented  model  of  the  missile  six-degree-of-freedom  dy¬ 
namics  describes  in  detail  the  motion  of  the  body  and  seeker  based  on  the  missile-target 
relative  kinematics.  However,  the  documented  model  is  not  based  on  the  standard  NED 
(North-East-Down)  axes  system  used  in  most  current  simulation  models  (Ref.  [1])  and,  in 
particular,  in  the  DREV  library  of  weapon  system  model  components.  In  order  to  main¬ 
tain  compatibility  with  the  existing  library  of  components,  the  documented  model  of  the 
missile  dynamics  has  been  modified  to  adopt  the  NED  standard  axes  system.  The  original 
and  NED  axes  systems  are  summarized  in  Table  I.  The  Euler  transformation  is  described 
in  Section  2.0.  The  modified  model  of  the  missile  dynamics  is  presented  in  detail  in  the 
following  sections. 


TABLE  I 


Definition  of  axes  systems 


Source 

Axes  definition 

Euler  transformations 

Original 

x  forward,  y  up,  z  lateral 

y,  zn,  Xi2 

NED 

x  forward,  y  lateral,  z  down 

z,yn,Xi2 

This  work  was  carried  out  under  Work  Unit  3ea20,  ”Air-to-Air  Missile  Model  De¬ 
velopment”  of  Project  3e,  ”Air  Weapon  Systems”,  from  March  to  May  1997. 
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In  the  standard  axes  system,  an  Euler  transformation  consists  of  three  consecutive  Euler 
rotations  around  the  z,  y  and  x  axes  respectively.  Each  successive  rotation  takes  place 
relative  to  the  intermediate  coordinate  system  resulting  from  the  previous  rotation.  The 
three  rotations  are  expressed  as  three  matrix  operators,  R$A[  ],  ReE[  ]  and  Ror\  ]>  which 
are  applied  in  sequence.  The  rotation  operator  for  6  a  represents  a  rotation  about  the  z 
axis,  and  is  defined  as 

cos(x)  sin(x)  0 

ReA[x]  =  —  sin(x)  cos(x)  0  [1] 

[  0  0  1 

The  rotation  matrix  for  6e  represents  a  rotation  about  the  yn  axis,  and  is  defined  as 

f  cos(x)  0  —  sin(x)  1 


ReE[x)  =010  12 

sin(x)  0  cos(x)  j 

and  the  rotation  matrix  for  6r  represents  a  rotation  about  the  X{2  axis,  and  is  defined  as 


ReR[A  = 


10  0 
0  cos(x)  sin(x) 
0  —  sin(x)  cos(x) 


The  global  transformation  operator  from  the  earth-fixed  coordinate  system  (x,  y,  z) 
to  any  coordinate  system  i  (xj,  j/i,  z»)  is  the  product  of  these  three  rotation  operators,  namely, 


Vi  =  ReR[d3,]  •  Roe[02\  ‘  V 


where  0l5  02  and  63  denote  the  three  successive  Euler  rotations.  The  inverse  transformation, 
from  coordinate  system  i  back  to  the  earth-fixed  system,  is 


y  =  Rsa[-6i]  ■  Rqb[-62]  ■  Ror[-03]  fji 


The  inverse  rotation  matrix  operators  are  denoted 

Roa[~Q l]  => 

Roe{-9 2]  =>  R-eB $2] 
RoR[-e 3]  => 


[6] 
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Using  the  missile  body  Euler  angles  referenced  to  the  earth-fixed  coordinate  system 
(namely,  6a,6e,6r)  as  an  example,  the  transformation  from  the  earth-fixed  to  the  mis¬ 
sile  body-fixed  system  consists  of  a  rotation  6az  followed  by  OeVh  and  finally  Qrx The 
coordinate  system  ( x ,  y,  z)  represents  the  initial  coordinate  system  for  the  transformation, 
(Xu, yn, Zi\)  and  (x^, y^, Z&)  denote  coordinate  systems  resulting  from  intermediate  Euler 
rotations,  6a  and  Be  respectively,  and  ( Xb,yb ,  Zb)  denotes  the  final  coordinate  system  result¬ 
ing  from  the  last  rotation,  Or.  Since  this  last  rotation  is  around  the  x&  axis,  x&  corresponds 
to  Xb ■  From  eq.  4,  the  global  transformation  in  this  example  is  summarized  as 


Xb 

X 

Vb 

.  . 

=  RoM-ReB[0E\-ReA[W 

y 

z 

[7] 


2.2  Coordinate  systems 


The  modeling  of  the  engagement  between  a  missile  and  its  target  consists  in 
describing  the  missile  and  target  dynamics  in  their  respective  coordinate  systems  and  their 
relative  kinematics  in  the  inertial  system.  A  minimum  of  six  different  coordinate  systems  are 
required,  and  they  are  summarized  in  Table  II.  The  earth-fixed  coordinate  system  represents 

TABLE  II 


Primary  coordinate  systems 


Type 

Acronym 

Axes 

Earth-fixed 

EF 

(z,  y,  z) 

Missile  body-fixed 

BF 

( Xb,Vb,Zb ) 

Missile  seeker-fixed 

SF 

(xs,ys,zs) 

Missile  velocity  vector-fixed 

VF 

{Xvi  Vv >  Zv) 

Missile-target  line-of-sight  fixed 

LF 

( Xi,yi,zi ) 

Target  body-fixed 

TF 

( *uVt,zt ) 

the  inertial  reference  frame.  The  missile  body-fixed  coordinate  system  is  attached  to  the 
missile  body,  and  its  x  axis  is  aligned  with  the  missile  longitudinal  axis,  the  y  axis  is  lateral 
and  the  z  axis  is  downwards.  The  missile  seeker-fixed  coordinate  system  is  positioned  at 
the  seeker  gyro  axis  of  rotation.  When  the  seeker  is  in  the  null  position  (boresight  along  the 
missile  body),  the  axes  of  the  seeker-fixed  system  are  aligned  parallel  to  those  of  the  missile 
body-fixed  system.  An  additional  coordinate  system  fixed  to  the  missile  velocity  vector  is 
required  to  compute  the  components  of  the  missile  incidence.  Also,  a  coordinate  system 
is  attached  to  the  missile-to-target  LOS  vector,  with  its  x  axis  along  the  LOS  vector,  to 
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compute  the  missile-target  relative  kinematics.  Finally,  a  coordinate  system  is  attached  to 
the  target  to  model  its  dynamics  and  kinematics. 

Euler  transformations  result  in  two  intermediate  coordinate  systems.  The  notation 
for  these  two  systems  is  summarized  in  Table  III.  For  simplicity,  this  notation  is  used  in  the 

TABLE  III 

Intermediate  coordinate  systems 


Type 

Notation 

Axes 

Intermediate  I 

11 

(®tlj  J/il)  Zil) 

Intermediate  II 

12 

(Xi2,Vi2,Zi2) 

present  document  to  identify  the  intermediate  coordinate  systems  of  all  Euler  transforma¬ 
tions.  Accordingly,  the  systems  (xu,yn,zn)  and  (xi2,  J/t2)  ^2)  refer  to  different  coordinate 
systems  for  different  Euler  transformations. 

The  Euler  transformations  between  the  various  coordinate  systems  of  Table  II  are 
summarized  in  Tables  IV,  V,  VI  and  VII,  where  the  Euler  rotations  are  shown  in  order. 
Table  IV  presents  the  Euler  transformations  from  the  earth-fixed  coordinate  system  to  each 
of  the  other  five  systems,  while  Thble  V  summarizes  the  transformations  from  the  missile 
body  to  the  seeker  and  missile  velocity  vector.  Table  VI  presents  the  transformation  from 
the  missile  body  intermediate  II  coordinate  system  to  the  missile  velocity  vector.  The  missile 
body  intermediate  II  system  is  obtained  by  performing  only  the  first  two  rotations  of  the 
inertial-to-missile  body  transformation  (eq.  7),  namely,  rotations  6a  and  Be  in  sequence. 
Finally,  Table  VII  presents  the  transformation  from  the  missile  seeker  to  the  missile-target 
LOS,  which  is  used  in  the  computation  of  the  seeker  tracking  error. 


TABLE  IV 

Euler  transformations  referenced  to  the  inertial  system 


Name 

Rotations 

missile  body 
seeker 

missile  velocity  vector 
LOS 

target  velocity  vector 

0AZ,dEyii,6RXi2  (xi2  is  xb) 
<f>AZ,<t>Eyn,<j>RXi2  (Xi2  is  xs) 
rYAZ,'yEyn,ryRXi2  (xi2  is  xv) 
0AZ,CTEyil,(7RXi2  {Xi2  is  &l ) 
TAZ,TEyiuTRXi2  ( Xi2  is  it) 
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TABLE  V 


Euler  transformations  referenced  to  the  missile  body-fixed  system 


Name 

Rotations 

missile  seeker 
missile  velocity  vector 

tRLZb,euDVn,  0xi2(xi2  is  x„) 

— OtRL&b ,  —OcuDVil,  0fj2  (Xi2  is  Xv) 

TABLE  VI 

Euler  transformations  referenced  to  the  missile  body  intermediate  II  system 


Name 

Rotations 

missile  velocity  vector 

-a A' Zb,  -(XEVii,  0x^2  (x{2  is  xv) 

Note:  The  intermediate  II  system  is  obtained  from  the  transformation  Rqe  [6e] -R$a  [9a\  ■  [£ , 


TABLE  VII 

Euler  transformations  referenced  to  the  missile  seeker  system 


Name 

Rotations 

LOS 

^RLZs,^UDVil  ,  0  Xi2  ( Xi2  is  x{) 
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The  Euler  transformations  of  Tables  IV,  V,  VI  and  VII  have  the  following  particu¬ 
larities: 

•  No  angular  £&  (xa)  rotation  from  the  body-fixed  to  the  seeker-fixed  systems  (Table 

V) 

•  No  angular  x &  (xv)  rotation  from  the  body-fixed  to  the  missile  velocity-fixed  systems 
(Table  V) 

•  No  angular  x &  (xv)  rotation  from  the  body-fixed  intermediate  II  to  the  missile 
velocity-fixed  systems  (Table  VI) 

•  No  angular  £&  (xi)  rotation  from  the  seeker-fixed  to  the  LOS-fixed  systems  (Table 
VII) 


The  non-Eulerian  angular  rates  of  the  missile  body,  missile  seeker  and  missile-target 
LOS  are  defined  in  Table  VIII.  These  rates  are  not  Euler  angular  rates  because  the  rotations 
take  place  around  the  three  axes  of  the  same  coordinate  system,  namely,  the  body-fixed, 
seeker-fixed  and  LOS-fixed  system  respectively. 


TABLE  VIII 

Non-Eulerian  angular  rates 


Name 

Angular  rates 

missile  body 

Orib^UDVb^RLZb 

missile  seeker 

<t>risAuDVsARLZs 

LOS 

&TXhOUDyi,ORLZl 

The  mathematical  relations  of  the  six-degree-of-freedom  missile  dynamics  model  are 
next  developed  using  the  standard  NED  axes  system  (Ref.  [1]  and  Table  I)  .  The  discussion 
also  includes  two  issues  of  importance  in  missile-target  engagement  modeling,  namely,  the 
computation  of  the  initial  conditions  and  of  the  miss  distance. 
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3.0  MISSILE  DYNAMICS 


3.1  Body  coordinate  system 

When  referenced  to  the  inertial  coordinate  system,  the  missile  body  angular  rate 
vector,  denoted  cD*,,  can  be  defined  in  terms  of  the  Eulerian  angular  rotations  9  a,  6e  and 
Or,  namely, 

tib  =  0RXi  2  +  0EVi  1  +  9az  [8] 

The  inertial  and  missile  body-fixed  coordinate  systems  are  related  by  the  Euler  transforma¬ 
tion  based  on  the  angles  0a,0e,0r  (eq.  7).  Equation  8  can  therefore  be  converted  to  the 
body-fixed  system  using  this  same  transformation.  The  transformation  from  the  (x,y,z) 
system  to  the  body-fixed  one  requires  three  rotations,.  Two  rotations  are  required  from 
(•Eii?  t/n,  zn)  to  (xb,yb,  Zb)  and  one  rotation  from  (x& ,  y& ,  2i2 )  •  Accordingly,  the  z  compo¬ 
nent  of  uJfc  in  eq.  8,  9 A,  must  be  rotated  through  all  three  rotations.  The  yn  component,  9e, 
must  be  rotated  through  two  rotations,  namely,  $e  and  6r.  The  £&  component,  6r,  only 
needs  to  be  rotated  through  Or.  In  summary,  the  transformation  of  the  three  components 
of  eq.  8  to  the  body-fixed  system  is 


0 

r  p  i 

Or 

^6  =  ReR  [9r\.R8e  [9e]-Rba  [9 a]- 

0 

6a 

+  Rbr[0r\-R6e[0e\- 

_ i 

+  R0r[Or]. 

0 

0 

Substituting  for  the  definition  of  the  transformation  operators  (eqs.  1,  2,  3)  in  eq.  9, 
ub  =  {-  sin  [9e]0a  +  0R}xb  +  {cos[<?e]  sin^j^A  +  cos[0fl]O£}y6+ 


{cos[0£;]  cosf^j^A  -  sin^j^Jift  [10] 

The  missile  body  angular  rate  vector  u>b  is  also  defined  in  the  body-fixed  coordinate  system 
as 

UJb  ~  0rXb  +  OuDVb  +  ORLZb  [11] 

Equating  the  Xb,yb,Zb  components  of  eqs.  10  and  11, 

0  A  —  sec[0£;](cos[0H]0itL  +  sin[0rt]0t/£))  [12] 

0e  =  —  sin[0/j]0jt£,  +  cos  [0r]0ud  [13] 

Or  =  syd[0e\0a  +  Or  [14] 


Equations  12,  13  and  14  define  the  missile  body  Euler  angular  rates  ( Oa,Oe,Or ) 
referenced  to  the  earth-fixed  coordinate  system  in  terms  of  the  non-Eulerian  rates 
( Or,OuD,0RL )  referenced  to  the  body-fixed  system.  The  body-fixed  rates  ( 0t,0ud,0rl ) 
are  computed  from  the  traditional  missile  equations  of  motion  (Section  3.5). 
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3.2  Seeker  coordinate  system 


The  seeker  angular  rate  vector,  denoted  ua,  is  expressed  in  terms  of  the  Euler 
angular  rotations  {(!>a,4>Ei4>e)  in  reference  to  the  earth-fixed  system  as 

W,  =  (j>RXi2  +  4>EVil  +  4>AZ  [15] 

us  is  also  expressed  in  the  seeker-fixed  coordinate  system  as 

UJ3  —  4>tXs  +  <j>UDVs  +  4>RLZs  [16] 


The  conversion  of  eq.  15  to  the  seeker-fixed  system  is  obtained  using  a  procedure  identical 
to  the  one  presented  in  Section  3.1  for  except  that  the  Euler  angles  are  now  (ft a,  <Pe  and 
4>r- 


=  Ror  [(fiR)-ReE[<fiE\- ReA  [4>a]  ■ 

0 

0 

+  ReR[<f>R]-ReE  I&e}- 

+  RorI^r]- 

<Pr 

0 

or, 

.  . 

L  0  J 

0 

[1 

tis  =  {4>r  ~  sm[4>E]4>A}xs  +  {cos [<f>E]  sin[(^ j?] <Pa  +  cos[<fiR]<fiE}ys 
-(-{cos[(^£']  cos[^]^x  -  sin [<Pr]4>e}zs 
Equating  the  Xg,ys,zs  components  of  eqs.  16  and  18, 


<f>A  =  secf^s]  (cos  [(Pr]4>rl  +  sin[<fiR]<fiUD)  [19] 

4>e~—  sin  [(I)r]4>rl  +  cos[(p  r]4>u  D  [20] 

cj)R  —  sin  [<Pe\^a  +  4>r  [21] 

Equations  19,  20  and  21  define  the  Euler  angular  rates  (4>Ai  4>Ei  4>r)  referenced  to  the  earth- 
fixed  system  in  terms  of  the  non-Eulerian  rates  ( <fir ,  <fiuDi  4>rl)  referenced  to  the  seeker-fixed 
system.  The  angles  <fiED  and  4>rl  are  computed  as  part  of  the  tracking  loop  (Section  4.0), 
while  <fir,  the  seeker  spin  rate,  is  obtained  in  Section  3.3. 


3.3  Seeker-to-body  coordinate  transformations 

This  Euler  transformation  is  more  complex  than  those  previously  discussed  because 
the  initial  coordinate  system  for  the  transformation,  namely,  the  body-fixed  coordinate 
system,  is  a  rotating  one  while,  previously,  the  earth-fixed  system  was  stationary.  The 
body  angular  rate  3b  is  related  to  the  seeker  angular  rate  us  according  to 


Ug  =  Cijb+Ua/b 


[22] 
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where  the  angular  rate  of  the  seeker  relative  to  the  body  is  denoted  as  uis/b-  The  terms  of 
eq.  22  are 

cSb  =  0rXb  +  OuDVb  +  &RLZb  [23] 

ujs  =  <prxs  +  (f>UDVs  +  4>RLZs  [24] 

and 

&s/b  =  ^UDVil  +  ^RLZb  [25] 

It  is  noted  that  there  are  no  x^  angular  rotation  from  the  body-fixed  coordinate  system  to 
the  seeker-fixed  system.  The  rotations  are  XfiLh  followed  by  XudVh  in  sequence.  Equations 
23,  24  and  25  are  converted  to  a  single  coordinate  system,  namely,  the  seeker-fixed  system. 
The  terms  in  eq.  23  refer  to  the  body-fixed  system  and  are  therefore  converted  to  the 
seeker-fixed  one  according  to 


£b  =  Rqe  [Xud\-Rqa  [Aitj,]. 


Or 

Qud 

Qrl 


[26] 


or 

{cos[AitL]  copfAt/D^j-  —  suifA^]#/^  +  cos[Aj/£>]  sin[A^]^i7£>}^«+ 

=  {—  sin[Aiti]0r  +  cos[A^£/]0f/i5}yJ+  [27] 

{cos[Afix]  sin[Ac/£>]^r  -I-  cos[Ai/c]0Ri,  +  sin[A/tL]  sin[A ud]Qud}zs 

The  term  us  requires  no  transformation.  The  term  ujs/b  is  transformed  to  the  seeker-fixed 
coordinate  system  using 


0 

r  o  i 

£s/b  =  R8e[Xud]-R6a[Xrl]- 

0 

.  ^RL 

+  Re  e[Xud\- 

— i 

Q 

fc>  o 
■-< 

[28] 


or 


w*/&  =  —  sin[A£/jo]AitLXs  +  A  udVs  +  cos[At/£>]Ajt£,£s  [29] 

Substituting  eqs.  24,  27  and  29  in  eq.  22,  and  equating  the  three  components  xs,  ys  and  zs, 


A  ud  —  sin[Aiti,]0r  —  cos[Aai,]0f/£)  +  <j>uD  [30] 

Xrl  =  -  sec[A*7£>](cos[Aitr,]  sin[A[/D]0r  +  cos[\UD]dRL+ 

sinfAiii,]  sin[A  udWud  —  <Prl) 

4>r  =  cos[Aitr,]  cos[Aj7£j]0r  —  sin[Aj/£>]0ft£  +  cosfA^]  sin[A  rl]0ud  —  sin[A[/£>]A.Rz,  [32] 

Xud,Xrl  and  <fir  describe  the  seeker  dynamics  relative  to  the  missile  body. 
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3.4  Velocity-fixed  coordinate  system 

The  x  axis  of  the  coordinate  system  fixed  to  the  missile  velocity  vector  Vm  is 
aligned  with  the  velocity  vector.  Also,  the  missile  velocity  vector  is  described  by  the  Euler 
angles  -ja,  1e  and  7#  referenced  to  the  earth-fixed  system.  Accordingly,  the  velocity  vector 
in  the  velocity-fixed  system  is  converted  to  the  earth-fixed  system  using  the  inverse  Euler 
rotations  —'Yr,  —7 e  and  —7 a  in  sequence.  Using  the  inverse  transformation  operators 
described  in  Section  2.1, 

Vm 

Vm  ~  R-eA['YA]-R-eB[lE\-R-eR['YR\-  0  [33] 

0 

or 

Vm  =  cos[7^]  cos[j e]Vmx  +  cosfrs]  sin^lV^y  -  ain[jE]VMz  [34] 

The  missile  body  Euler  angles  referenced  to  the  earth-fixed  system  axe  (0a,  6e,  Or). 
Rotations  through  6a  and  Or  in  sequence  result  in  the  intermediate  reference  system 
{xi2,yi2,Zi2),  where  x&  corresponds  to  the  missile  body  longitudinal  axis,  The  an- 
gle  between  the  missile  velocity  vector  Vm  and  the  body  projected  on  the  x& ,  ya  plane  is 
defined  as  a  a1-  The  angle  perpendicular  to  the  x^ ,  m2  plane  is  defined  asag.  a  a  is  defined 
as  the  projection  of  a  a1  onto  the  earth-fixed  x,  y  plane.  A  final  rotation  Or  about  x&  results 
in  the  body  coordinate  system  (xb,yb,Zb).  aRL  is  defined  as  the  angle  between  the  velocity 
vector  and  the  missile  body  in  the  Xb,yb  plane,  and  auD  is  the  angle  perpendicular  to  the 
Xb,Vb  plane. 

Using  the  above  definitions,  the  missile  velocity  vector  can  be  expressed  in  both  the 
(xb,  yb,  Zb)  and  (xi2,  y&,  z&)  coordinate  systems.  Referenced  to  the  body,  the  velocity  vector 
is 

Vm 

Vm  =  R-eA[-(XRL]-R-6E[-auD]-  0  [35] 

0 

The  transformation  to  (xi2,y&,  ^2)  from  the  earth-fixed  system  is  performed  using  the  op¬ 
erator  R0e[Oe]-R0a[Oa],  and  the  last  rotation,  Rsr[0r],  converts  to  ( Xb,yb,Zb )■  Accordingly, 
the  complete  transformation  of  the  velocity  vector  through  (xb,yb,Zb)  to  is 

'  VM  ' 

Vm  —  R-6R[0R\-R-BA\-<XRL]-R-6E[~aUD]-  0  [36] 

0 
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°r 

VM  ~  cos[aK£,]  cos[aUD]VMXi2+ 

(-  cosfat/c]  cos[0R]  sin[a.Ri,]  -  sinfatrp]  sin[6R])VMyi2+  [37] 

(cos[0fl]  sin[a[7£>]  -  cos[auD]  sin[a#/,]  sin[0jt])VA/;Zt2 

It  is  noted  that  the  negative  of  the  incidence  angles  am,  and  ajjD  are  used  in  the  above 
transformation  since  the  incidence  is  measured  from  the  vector  Vm  to  the  missile  body.  The 
inverse  rotation,  from  the  body  to  Vm,  therefore  uses  negative  a' s.  Similarly,  the  velocity 
vector  referenced  to  the  system  {xi2,Vi2,  za)  is 


Vm  =  R-eA[-OLA']-R-eE[-<XE]- 


VM 

0 

0 


or 

VM  -  cos[q!£;]  cos[aA']VMXi2  -  cosfcnjs]  sm[aA>]VMyi2  +  sm[aE]VMZi2 

The  angles  aA>  and  aR  are  solved  from  eqs.  37  and  39, 


[38] 


[39] 


sin[a£;]  =  cos[0fl]  sinfa^r)]  -  cos[aj7/)]  sinfa/tL]  sin[0^] 


and 


sin[d!A']  =  sec[aE](cos[a^£)]  cos[dR]  sin[aKL]  +  sinfaf/^]  sin[0R]) 
From  geometrical  considerations, 

sinforyj']  =  sinfo^]  cos[7e] 


[40] 

[41] 

[42] 


If  the  missile  incidence  components  auD,  ocm,,  oe  and  aA>  remain  small  (less  than 
10  degrees),  eqs.  40  and  41  can  be  simplified,  after  incorporating  eq.  42,  to 

aR  —  —  sin[0/e]aiii  +  cos  [0R]auD  [43] 

and 

aA  =  sec[7E](cos[0fi]o!itL  +  sin^auo)  [44] 

where  jA,  je  axe  the  velocity  vector  Euler  angles  relative  to  the  earth-fixed  system, 

7a  —  9  A  ~  a  A 


7 e  =  0e~  cue 


[45] 

[46] 
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3.5  Equations  of  motion 

3.5.1  Translational  equations 


The  equations  of  motions  in  translation  are 

Fx  =  Mm  (u  —  Orlv  +  Oudw) 
Fy  =  Mm( v  +  6rlu  -  6rw) 
Fz  =  Mm{w  +  OrV  —  Ojjdv) 


Also,  since 
then, 

The  acceleration  is  given  by 


VM  —  uxb  +  vyb  +  wzb 


VM  =  \Vm 


Vm  = 

Solving  for  u,v,  w  from  eqs.  47,  48,  49, 


=  Vu2  +  v2  +  w2 

(uii  4-  vv  +  ww) 

V^i 


u  = 


Fx  +  MmvQrl  —  Mmw6ud 
Mm 


Fy  +  Mmw9t  —  MmuQrl 
V~  Mm 

Fz  —  Mmv9t  +  MmuGud 
w  = - 


Mm 


Substituting  eqs.  53,  54  and  55  into  eq.  52, 


■  _  Fxu  +  Fy v  +  Fz w 

M  ~  MmVm 


[47] 

[48] 

[49] 

[50] 

[51] 

[52] 


[53] 

[54] 

[55] 


[56] 


Since  the  rotation  from  the  body  to  the  missile  velocity  vector  Vm  is  —curl  followed 
by  —auD,  the  velocity  vector  is  converted  to  the  body-fixed  system  using  the  transformation 
Rbe[— OiiL]  •  ReA[—otuD ],  or 

Vm  =  cosfaftzj  cos[auD]VMXb  —  cos[aj/£>]  sin[aat,]  VMVb  +  sm[auD]VMZb  [57] 
Using  small  angle  approximations, 

Vm  =  VM%b  -  VMOCRLVb  +  VMOCUDh 


[58] 
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The  values  of  u,  u,  w  are  obtained  by  comparing  eqs.  50  and  58.  Substituting  in  eq.  56, 


T>  Fx  ~  FyOLRL  +  Fzoljjd 

Vm  = - HZ - 


The  components  of  the  missile  incidence  are 


CtRL  =  ~ 


VM 


and 


auD  = 


w 

Vm 


[59] 

[60] 
[61] 


Equations  60  and  61  are  differentiated  with  respect  to  time  to  obtain  expressions  for  am 
and  auD-  Substituting  eq.  54  for  v,  eq.  55  for  w  and  eq.  56  for  Vm, 

Fv  FxaRL 


and 


&RL  =  ORL  ~ 

ocud  =  Oud  + 


MmVm  MmVm 
MmVm  +  ^rQ 


; —  0rauD 

f 

[62] 

FxauD 

MmVm 

[63] 

3.5.2  Rotational  equations 

Given  that  the  missile  is  symmetric  and  neglecting  all  products  of  inertia  terms, 


Or  —  VI £  / 1 xx 

[64] 

Oud  =  {My  {I XX  I zz)0  RlOT)  /  Iyy 

[65] 

OrJL,  =  {Mz  —  {Iyy  Ixx)0ud0t) / Izz 

[66] 

where  Mx,My,Mz  denote  the  aerodynamic  moments  and  Ixx,  Iyy,  Izz, 
moments  of  inertia  about  the  three  primary  body  axes. 

the  missile  airframe 

3.6  Accelerometer  output 


The  conversion  of  the  gravity  force  g  from  the  earth-fixed  reference  system  to  the 
body-fixed  one  is 


9  =  ReR  [Or] -ReE  [9e\-Rba  [8  a]  •  {0, 0,g} 


[67] 
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or 


g  =  -g  sin^jfi  +  g  cos[0£]  sin[0fl]y6  +  g  cos[0s]  cos[0h]z6  [68] 


The  components  of  the  gravity  force  in  the  body-fixed  system  axe  subtracted  from 
the  net  missile  body  force  components  to  reflect  the  output  of  the  accelerometer.  The 
longitudinal  and  lateral  accelerations  measured  by  the  accelerometer,  in  units  of  g’s,  are 


r 

r  Fx+gsm\6E]MM  -I 

axb 

9mm 

Fv-9cosf0Elsin[0R].WAf  di„6nr. 

aVb 

.  a*»  . 

g. 

Fz-gcos]eE]cos]eR]MM  .  di„0nD 

L  jMm  9  -1 

where  dia  is  the  distance  from  the  center  of  gravity  to  the  accelerometer. 
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4.0  TRACKING  LOOP  DYNAMICS 


4.1  Target  kinematics 


The  target-velocity  Vt  is  assumed  to  be  aligned  with  its  body  at  all  times  during 
its  flight.  The  x  axis  of  the  target-fixed  coordinate  system  is  therefore  aligned  with  Vt- 
The  target  velocity  is  also  described  by  the  Euler  angles  U4,  T#  and  Ur,  referenced  to  the 
earth-fixed  system.  Accordingly,  the  target  velocity  in  the  target-fixed  system  is  converted 
to  the  earth-fixed  system  using  the  inverse  Euler  rotations  —  Ur,  —Te  and  —  Ta  in  sequence. 
Using  the  inverse  transformation  operators  described  in  Section  2.1, 


or 


VT 

0 

0 


[70] 


Vt  =  cosjU/i]  cos[T£;]Ur£  +  cosp^] sinjU^Vry  -  sin[r£]Vrz  [71] 


4.2  LOS  Euler  angles  and  angular  rates 
4.2.1  LOS  Euler  angles 

The  Euler  angles  of  the  missile- target  LOS  vector  are  computed  from  the  missile- 
target  relative  geometry  in  the  earth-fixed  coordinate  system.  The  relative  missile- 
target  range  vector  Rmt  has  components  (%mt,  yMT,  zmt)  in  the  inertial  system  and 
{xmt-IiVmt-i-,  zmt-i)  in  the  LOS-fixed  system.  Since  the  xi  axis  is  aligned  with  the  LOS 
vector,  the  yi  and  zi  components  of  the  range  vector,  namely  Vmt-i,  zmt-1,  are  null  and 
the  xi  component  is  the  missile-target  range,  Rmt-  Accordingly,  the  range  vector  is,  by 
inspection, 

Rmt  =  cos  [cr^]  cos[a  e]Rmtx  +  cos[cte]  sin[c a\RmtV  —  s'm[&E\RMTZ  —  Rmt&i  [72] 
The  Euler  angles  a  a-,  &e  are  obtained  by  inspection  of  the  geometry, 

OA~  tan- 1  [———■]  [73] 

XMT 
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The  coordinates  of  the  LOS  vector  in  the  LOS-fixed  coordinate  system  are  obtained  from 


xMT—l 

VMT-l 

=  Rqr  kfl] -Rbe^e)-  R#a  [*a]  ■ 

XMT 

Vmt 

_  xMT-l  _ 

.  ZMT 

[75] 


where,  for  instance,  is  xt~~xm-  Substituting  eqs.  73  and  74  for  a  a  and  aE  respectively 
in  eq.  75, 

xMT—i  =  Rmt,  Vmt-i  —  0,  zmt-i  =  0  [76] 


as  expected. 


The  expression  for  aE  is  obtained  from  the  earth-fixed  to  LOS-fixed  Euler  transfor¬ 
mation.  This  transformation  is  written  as  a  combination  of  two  transformations,  from  the 
earth-fixed  system  to  the  seeker-fixed  one  followed  by  one  from  the  seeker-fixed  system  to 
the  LOS-fixed  one,  namely, 

T\  =  RgE[euD]-R6A[eRL]-RoR[<f>R]-R6Ey>E\-R6A[<f>A]  [77] 

The  transformation  from  the  earth-fixed  to  the  LOS-fixed  coordinate  systems  can  also  be 
written  in  terms  of  the  LOS  Euler  angles  referenced  to  the  earth-fixed  system,  namely, 

Tb  =  ReR[<xR]-R0E[<7E]-R8AWA]  [78] 

From  eq.  78,  the  ratio  of  T2(2, 3)  over  T2( 3,3)  gives  tan[(jR]  since  the  two  terms  are 
cos[cte]  sin[cr.R]  and  cos [cr^]  cos[<xr]  respectively  (derivation  not  shown).  Since  Ti  (eq.  77) 
and  T2  (eq.  78)  are  equivalent,  the  complete  expression  for  tan[er/j]  is  given  by  the  ratio  of 
Ti  (2, 3)  over  Tx  (3, 3),  namely, 

tan[cr*]  =  (sin[e/2i]  sin^s]  +  cos[e  rl]  cos[(f)E]  sin[<fo])/ 

(-  cos[efl£]  sin[et/xj]  sin[<£B]  +  cosf^s] 

(cos[ecr£»]  cos[<f>R]  +  sinfeat,]  sin[et/£>]  sin[<f>R])) 

Equation  79  can  be  simplified  for  small  angles  tm,  and  euD, 

1  _  cos{<j)E]  $m[<pR]  +  sm[4>E]eRL 
an  CTH  cos[^js]  cos[<£;j]  -  8m[cf)E]euD 


[79] 


[80] 


4.2.2  LOS  Eulerian  angular  rates 


The  relative  missile-target  velocity  vector  Vmt  in  the  earth-fixed  system  is  defined 


as 


Vmt  —  xmtx  +  VmtV  +  zmtz 


[81] 
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The  relative  velocity  vector  referenced  to  the  earth-fixed  system  can  be  converted  to  the 
LOS-fixed  system  using 


Vmt  =  Ror  [°b\-Rge  We]  Re a  Wa] • 


&MT 

Vmt 

zmt 


[82] 


or 

Vmt  =  {cos^j  cosWe]&mt  +  cos[<te]  s\n[a A]yMT  ~  8m[aE]zMT}xi+ 

{(-  cos  [cr^j]  sin[o-A]  +  cos[<ta]  sin[<7£]  sin[cr/j])iMT  +  (cos[<ta]  cos[<xh]  +  sin[<7A]  sin[er£] 
sin[<rK])yA/r  +  cos[crE]  sin  [aR]zMT}yi+ 

{(cos[cta]  cos[ctr]  sinfcr^]  +  sin[ovi]  sin[aR])xMT  +  (cos[ctr]  sinfcr^]  sinfo-jj]  —  cos[<xa] 
sin[CTfi])yMr  +  cosfcrjs]  cos[aR]zMT}Zi 

[83] 


The  missile-target  velocity  vector  relative  to  the  earth-fixed  system  is  also  obtained 
by  differentiating  eq.  72  with  respect  to  time  t, 

Vmt  =  Rmt  =  {cos[<ta]  cos[ge]RMt  -  cos[ct£]  sin Wa]Rmt*a  ~  cos[cta]  sm[aE]RMTcrE}x+ 
{cos[<7£]  s'm[<7 a}RMt  +  cos  [cta]  cos{ge]Rmt^a  -  sin[cM]  sm[aE]RMT^E}y+ 
{—  s\h\ge]Rmt  —  cos{oe}Rmt&e}z 

[84] 

Equation  84  is  converted  to  the  LOS-fixed  coordinate  system  using  the  Euler  transformation 

based  on  the  three  Euler  rotations  aA ,  a E  and  ctr,  and  the  result  is 

— *  * 

Vmt  =  Rmt^i 

+Rmt{cos[oe]  cos[cth]cta  -  sm[aR]dE)yi  [85] 

--Rmt(cos[£T£]  sin[cr/d]crA  +  COsWrWe)zi 

where,  as  expected,  the  xi  component  is  Rmt •  Since  eqs.  83  and  85  are  equivalent,  their 
( &i,yi ,  zi)  components  are  equated.  The  xi  term  is  denoted  the  closing  velocity  and  is  given 
by 

Rmt  -  cos[cta]  cosWeWmt  +  cos[a£]  s\n[a A\yMT  ~  sm[aE\zMT  [86] 

The  yi  and  zi  components  are  solved  simultaneously  for  aA  and  aE,  namely, 

a  a  =  -  •— — {sec  [erg]  sm[aA]xMT  ~  cos[cta]  sec[a  E]yMT}  [87] 

Rmt 

aE  =  — rri— (cos[<ta]  sin[crf;]xMr  +  sin[<TA]  sin[cr£:]yMT  +  cos[ge]zmt)  [88] 
Rmt 

The  simultaneous  solution  of  the  yi  and  zi  components  removes  any  dependence  on  aR. 


The  third  component  of  the  LOS  angular  rate  vector,  ctr,  is  obtained  as  follows. 
The  angular  rate  vector  of  the  LOS,  u y,  can  be  defined  as 


<2i  =  ORXi2  4-  aExn  +  aAx 


[89] 
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It  is  converted  to  the  LOS-fixed  system  using 

0  0  OR 

-  ReR[°R]-RoB[oE]-ReA[cA]-  0  +i^*[crH]--R0s[crE]-  oE  +R6r[or\.  0  [90] 

a  a  0  0 


or 

u>l=  —  sin[crB]d/i  4-  orxi + 

cos[£7S]  sin[cxii]o'/l  +  cos[or]&eVi+  [91] 

cos [<7£]  cos^it]^  -  sm[oR]oEzi 

Cbi  is  also  expressed  in  terms  of  the  non-Eulerian  angular  rate  components  relative  to  the 
earth-fixed  reference  system  (wie,wfe,wie),  as 


+  ujyey  +  uzez  = 


Rmt  x  Vmt 

&MT 


Converting  eq.  92  to  the  LOS-fixed  reference  system, 


&XC 

til  —  R6r[o-r].R6b[(Te]-R6aWa]-  Wye  [93] 

_  Vze  _ 

Equating  £[  components  of  eqs.  91  and  93  and  substituting  for  a  a  from  eq.  87, 
or  =  cos[cr a]  cos[oE]ujxe+ 

cos[cr£]  sin [oA}uye  -  sin [oE]u)ze  -  sin[crE]  [94] 

(sec[crjs]  svd{oa]xmt  ~  cos  [a  a]  sec[oE]yMT) 


The  seeker  non-Eulerian  rotation  rates  in  the  seeker-fixed  system,  or ,  oud  and  orl, 
are  determined  in  terms  of  the  Eulerian  rates  oa,  oe,  or,  using  the  approach  presented  in 
Section  3.1.  The  Eulerian  angular  rotation  vector  referenced  to  the  earth-fixed  system  is 
given  by 

Cji  =  ORXi  2  +  OEVn  +  °AZ  [95] 

The  angular  rotation  vector  expressed  in  the  LOS-fixed  system  is 

Ui  =  orxi  +  ouDm  +  orlzi  [96] 

Equation  95  is  converted  to  the  LOS-fixed  coordinate  system: 

0  0  or 

3i  =  RoR[0R\-ReB[aE]>RoA[oA)'  0  +Ror[or].Rob  [oe\-  oe  +Rbr[or ].  0  [97] 

&A  0  0 
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or 

U \  =  (-  sin[crB](TA  +  crR)xi  +  (cos[cth]  sin [or\ca  +  cos[cr  r\ct  E)yi+ 

[98] 

(cos[ct£;]  cos  \orJ<j  a  —  am[api\dE)zi 

Equating  the  yi,  zi  components  of  eqs.  96  and  98, 

CTRL  —  cos [cth]  COS [ctr] <7 a  —  $in[a r\oe 

[99] 

ctvd  —  cos[cth]  &m[<jR\pA  +  cos[u  r}ge 

[ioo] 

or  —  —  siu[oe]oa  +  OR 

[101] 

The  above  relations  describe  the  LOS  angular  rates  or,&RL,ouD  referenced  to  the 
LOS-fixed  system  in  terms  of  the  LOS  Euler  angular  rates  &a,  or,  or  referenced  to  the 
earth-fixed  system. 

The  angular  separation  between  the  LOS  vector  and  the  gyro  spin  axis  (seeker 
boresight)  is  denoted  e  and  its  components  are  €rl  and  euD-  The  components  represent 
Euler  rotations  in  sequence  from  the  seeker-fixed  system  to  the  LOS-fixed  one,  namely, 

e  =  euDVil  +  tRLZs  [102] 

The  total  tracking  error,  e,  is  defined  as 

e  =  \jeRL  +  eUD  [103] 

The  tracking  rate  4>  is  proportional  to  the  total  tracking  error.  As  a  first  approxi¬ 
mation,  the  tracking  rate  (f>  can  be  modeled  in  a  generic  way  using 

4>  =  KpeH(s)  [104] 

where  H(s )  is  the  seeker  tracking  loop  transfer  function  and  Kp  is  a  constant  gain. 

— * 

The  direction  of  the  tracking  rate  vector,  4>,  is  given  by  xs  x  xi,  that  is,  along  the 
vector  perpendicular  to  the  plane  formed  by  xa  and  xi-  The  unit  vector  xp  is  defined  to  be 
in  this  direction,  namely, 

„  _  Xs  X  Xl 

Xp  \x3  X  Xl\ 

$  is  hence 

<j>  —  KpeH(s)xp  [106] 

The  term  xa  x  xi  is  obtained  from  the  seeker-fixed  to  the  LOS-fixed  Euler  transfor¬ 
mation  {ejjDitRL)-  The  xi  component  of  the  transformation  R$E  [ey d\-Roa  [crl]  is 


xi  =  cos [e Hi]  cos [e£/p]xs  4-  cos[ej7i>]  sin[eRi]j/.,  -  sin[e  £/£>]£,, 


[107] 
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Accordingly, 

xsx£[  =  sin  [euD]Vs  +  cos[epp]  sinfe/i/Ji, 

[108] 

and 

\x„  x  &i\2  =  cos [ci/d] 2sin[eitt] 2  4-  sin[et/£>]2 

[109] 

For  small  angles,  cos[e 

]  w  1,  and  sin[e]  «  e, 

x»  x  xi  fa  eUDys  +  c-RLZs 

[110] 

and,  using  eq.  103, 

|s4  x  s/|  w  e 

[111] 

Accordingly,  eq.  105  for  xp  reduces  to 

CUD  -  ,  £rl  . 

Xp  =  e  Vs  +  e  Zs 

[112] 

Finally,  substituting  eq.  112  in  eq.  106,  $  referenced  to  the  seeker-fixed  system  is 

$  =  KpH(s)  [ej jdVs  +  f-RLZa]  [113] 

The  components  of  (p,  denoted  <puD  and  4>rl,  are 

— ♦  •  • 

<t>  =  <t>UDVs  +  <pRLZs  [114] 

In  comparison,  the  seeker  non-Eulerian  angular  rate  referenced  to  the  seeker-fixed  system 
is 

us  =  4>rxs  +  pUDVs  +  4>rlzs  [115] 

The  rates  €rl  and  euD  are  derived  using  the  procedure  outlined  in  Section  3.3.  The 
LOS  angular  rate  vector  ui  is  defined  in  reference  to  the  seeker-fixed  system  as 

&i  =  <jjs  +  <2i/s  [116] 

The  LOS  non-Eulerian  angular  rate  referenced  to  the  LOS-fixed  system  is 

uji  =  arxi  +  oudVi  +  <?RLZt  [117] 

The  Euler ian  angular  rate  vector  of  the  LOS  referenced  to  the  seeker-fixed  system  is  denoted 
as  Cbi/g  and  is  defined  as 


&l/s  =  eUDVil  +  tRLZs 


[118] 
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Equations  115,  117  and  118  axe  converted  to  a  single  coordinate  system,  namely,  the  LOS- 
fixed  system.  ujs  is  converted  using  the  transformation  R$E [euD]-ReA [crl]i  namely, 

Cos  =  {cos[ejtt,j  cos[ei/£>]^r  ~^[(-ud]4>rl  +  cos[ej7£>]  sin [€rl]4>ud}£i+ 

{-  sin[eat,]^r  +  cos[e  RL]4>UD}yi+  _  [119] 

{cos[e«L]  sm[euD]<j>r  +cos  [cud^rl  +  sinfc^]  sm[eUD\4>UD}zi 


The  LOS  angular  rate  vector  ui  is  already  in  the  LF  system.  Coi/s  is  obtained  from  the 
transformation 


or 


£i/a  —  ReE{euD]-ReA[cRL]- 

0 

0 

+  R6e[€Ud] ■ 

0 

CUD 

crl 

0 

[120] 


&i/s  =  -  sin [euD^RLXi  +  kuDVl  +  cos [cud]crlzi 


[121] 


Equations  117,  119  and  121  are  substituted  in  eq.  116.  Equating  the  yi  components 
and,  given  that  the  maximum  angular  separation  e  is  usually  less  than  10  degrees,  the 
resulting  expression  for  kuD  is 

CUD  =  &UD  +  CRl4>t  —  4>ud  [122] 

Similarly,  equating  the  zi  components, 

C-RL  =  CTRL  -  euD<j*r  -  <f>RL  [123] 

Equations  122  and  123  are  integrated  to  give  crl  and  cud  as  a  function  of  time.  The  total 
error  e  is  computed  from  eq.  103  and  substituted  in  eq.  113  to  give  the  total  tracking  rate 
vector  <£  (specifically,  ud,4>rl )• 

The  total  tracking  rate  will  usually  be  limited,  and  the  limited  value  is  denoted  as 

• 

<f>i.  The  components  of  4>i  are 

$i  =  4>i%p  [124] 

7  eUD  \  -  ■  €RL x  ~ 

<f>l  =  —^-VlVs  +  —Wz* 


or,  using  eq.  112, 


[125] 
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5.0  SIMULATION  CONTROL 


5.1  Initial  conditions 

The  initial  conditions  are  obtained  from  the  following  specified  parameters: 


•  Vm,Vt 

•  °A,  E 

•  9  a  ,  0e  ,  Or 

•  rA,rE  (tr  =  o) 

•  Rmt 

•  hx  (target  altitude) 

•  1A,7E  {lR  =  0) 

The  initial  components  of  the  missile  incidence,  aA,aE  are  obtained  from  eqs.  45 
and  46.  The  relations  for  the  initial  values  of  oiud,<xrl,  ^ud,  ^RLi  4>Ai  <t>E  and  4>r  are 
summarized.  The  angles  olud  and  (xrl  are  computed  from  eqs.  43  and  44, 

aUD  =  cos[7£]  sin[^]aA  +  cos[0e]c*£  [126] 

and 

ctRL  =  cos[7£]  cos^/jjaA  -  sin^ttE  [127] 


Seeker  lock-on  is  assumed  at  missile  launch.  The  gyro  axis  xs  is  hence  aligned  with 
the  LOS  vector  and  thus  corresponds  to  x/.  xi  is  given  by  the  following  rotation  from  the 
earth-fixed  system, 

r  *,  i  i  *  i 

[128] 


Xl 

X 

Vl 

=  ReR  Wr] -r6e  [oe] lie  a  M 

y 

z 

The  transformation  from  the  earth-fixed  system  to  the  body-fixed  one  is  given  by 


Xb 

X 

Vb 

.  . 

=  r6r  [0ft]  •  ReE  [9e] -r6a  [9a] 

V 

z 

[129] 
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The  unit  vector  xt,  is  therefore  transformed  to  the  LOS-fixed  system  using 


Xl\ib 

xb 

Vl\ib 

—  ReR  [<tr].Rob  [oe\-Rqa  [<7a]-R-6a  [Qa]-R-oe  [9E]-R-eR  [0fi]- 

0 

0  _ 

[130] 


The  contribution  of  Xb  to  xi,  denoted  Xi|x6,  is,  from  eq.  130, 

xi  \xb  =  (cos  [0^]  (cos  [0a]  cos  [a  a]  cos[o\e]  +  cos  [ob]  8in[0A]  sin[<7A])+ 
sin[0B]  sin[<T£:])xfc 

Similarly,  the  contribution  of  yb  to  xi  (eq.  130),  denoted  as  is 


[131] 


Xi  |j/6  =  (cos  [0b]  (-  cos  [cr  a]  cos[ob]  sin[0A]  +  cos[0a]  cos[(7b]  sin[(7A])  +  sin[0n] 

(sin[0B](cos[0A]  cos [(7a]  cos[(7b]  +  cos [ctb]  sin[0A]  sin[(7A])  -  cos [0b]  sin[(7B]))y& 

[132] 

and  the  contribution  of  Zb  to  xi  (eq.  130),  denoted  as  xi\zb,  is 

xi\ib=  (-  sin  [0  b]  (-  cos  [(7a]  cos  [ctb]  sin[0A]  +  cos[0a]  cos  [ctb]  sin[(7A])  +  cos  [0b] 

(sin[0B](cos[0A]  cos [(7a]  cos[(7b]  +  cos[ob]  sin[0A]  sin[(7A])  -  cos[0b]  8in[<rB]))2& 

[133] 

The  contributions  of  Xb ,  yb,  %b  to  are  also  obtained  from  the  missile  body  to  seeker 
transformation.  The  transformation  of  Xb  to  the  seeker-fixed  system  is  performed  using  the 
transformation 

T  xs\xh  1  r  Xb  1 

[134] 

This  is  also  the  contribution  to  the  LOS-fixed  system  since  the  two  systems  are  initially 
co-located: 

r  l  r  xs\£k  i 

[135] 


Xs\xb 

Xb 

Vs\xb 

—  R6e[^ud]-R0a[^rl\- 

0 

Zs\tb 

0 

%l\xb 

Vlkb 

= 

Vs\xb 

Zl\xb 

1 

H> 

cr 

1 _ 

Accordingly,  the  contribution  of  Xb  to  xi  (eq.  134)  is 

=  COS  [Abb]  cos[A{/b]^6 

Similarly,  the  contribution  of  yb  to  xi  (eq.  134)  is 

xi\yb  =  cos[Abb]  sin[ABL]y& 
and  the  contribution  of  Zb  to  xi  (eq.  134)  is 

xi\ab  =  -sin[Aj7D]4 


[136] 


[137] 


[138] 
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From  the  eqs.  136,  137  and  138,  — equals  sinfAt/zj]  and  the  ratio  of  to  &i\ib  equals 
tan[Aj?L]. 

Since  eqs.  131,  132  and  133  correspond  to  eqs.  136,  137  and  138  respectively, 
sinfAt/zj]  is  given  by 

sin[A{/£)]  =  sin[0fl](—  cos  [ax]  cos[cr£:]  sin[#x]  +  cos[0x]  cos[<7£;]  sin[ax])  —  cos[#/j] 

(sin[0£;](cos[0x]  cos  [ax]  cos[as]  +  cos[ajg]  sin[#x]  sin[ax])  —  cos[0£]  sinfa#]) 

[139] 

and  tan[A^],  by 

tanfAnz,]  =  (cos  [$#](—  cos  [ax]  cos[az;]  sin[#x]  4-  cos[#x]  cosfcr^]  sin[ax])  4-  sinf##] 

(sin[0£](cos[#x]  cos[ax]  cos[a£?]  4-  cos[a£:]  sin[#x]  sin[ax])  —  cos[#£]  sin[a£]))/ 
(cos[0£’](cos[0J4]  cos  [ax]  cos[a#]  4-  cosfaz;]  sin[#x]  sinfax])  4-  sinf##]  sinfae]) 

[140] 

The  orientation  of  the  seeker  coordinate  system  relative  to  the  earth-fixed  system 
is  given  by  the  Euler  angles  </>x,  <Ae,  4>r-  Since  the  seeker-fixed  and  LOS-fixed  systems  are 
co-located  at  launch,  the  seeker  boresight  error  angles,  euD,  crl,  are  null  and  the  seeker 
boresight  Euler  angles  <^>x  and  4>e  are 


<t>A  =  <?A 


and 

relative  to  the  earth-fixed  system. 


4>e  —  <?E 


[141] 

[142] 


The  initial  roll  angle  <j>R  is  obtained  next  using  two  expressions,  one  derived  from 
the  LOS-fixed  system  and  a  second  one  from  the  seeker-fixed  system.  Seeker  lock-on  at 
launch  implies  that  <f>R,  the  Euler  roll  angle  of  the  seeker-fixed  system,  equals  or,  the  Euler 
roll  angle  of  the  LOS-fixed  system.  First,  the  intermediate  II  coordinate  system  of  the 
Euler  transformation  from  the  earth-fixed  system  to  the  LOS-fixed  one  is  considered.  The 
intermediate  II  system  (£j2,  Vi2,  za)  is  obtained  by  rotating  the  earth-fixed  system  (5,  y,  z) 
through  the  Euler  rotations  a  a  and  a#  in  sequence.  The  y*2  term  is 

Vi2  =  -  sin[ax]£  4-  cos[ax]y  [143] 

The  unit  vector  y» 2  is  now  expressed  in  the  body-fixed  reference  system  by  first  rotating 
back  to  the  earth-fixed  system  and  forward  to  the  body-fixed  system  according  to 


r  °  1 

Vb\yi2 

^lv<2  . 

=  Rbr  [fo] -ReE  [8e}-Roa  [0a]  R-ba  [va]-R-be  [&e}- 

- 1 

<N  _ 

<S° 

_ 

[144] 
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or 


{cos[0e]  cos[cta]  sin[$A]  —  cos  [0  a]  cos[0ft]  sin[crA]}> 

^b\ji2 

{cos[cta](cos[0a]  cos[0ft]  +  sin[0A]  sin[0£]  sin[0ft])- 

Vb\yi2 

= 

(-  cos[0ft]  sin[0A]  +  cos  [0a]  sin[0ft]  sin[0ft])  sin[tJA]}, 

*b\i/i2 

(cos[cta](cos[0/{]  sin[0A]  sin[0g]  -  cos[0a]  sin[0ft])- 

.  (cos[0a]  cos [0ft]  sin[0ft]  +  sin[0A]  sin[0ft])  sinfoA]} 

from  which  the  first  expression  for  ya-h  is  obtained. 


[145] 


A  second  expression  for  ya-Zb  is  obtained  by  considering  the  seeker-fixed  system. 
The  transformation  from  the  seeker-fixed  system  to  the  body-fixed  one  is 


xb 

Xs 

Vb 

.  *h  . 

—  R-oa[\rl].R-8b{\vd]- 

Vs 

.  *s  . 

where  the  Zb  term  is 

Zb  —  —  sin[At/p]xs  +  cosfAj/cjz, 


[146] 


[147] 


The  transformation  from  the  seeker-fixed  to  the  intermediate  II  system  is  only  i2_0fl[0ft], 
or  equivalently,  i?_0fi[<7ft],  as  noted  previously.  The  ya  component  is  hence 


yi2  =  cos [<f>R]ys  -  sin[^]i5  [148] 

The  second  expression  for  yi2.Zb  is  obtained  from  eqs.  147  and  148,  namely, 


Va-h  -  -  cos[Ac/d]  sin[<fo] 

Combining  with  eq.  145  and  solving  for  sin[</>ft], 

sin[0fl]  =  -  sec[At/p](cos[cr/i](cos[0fl]  sin  [0  a]  sin[0ft]  -  cos [0a]  sin[0ft])— 
(cos[0a]  cos[0ft]  sin[0£]  +  sin[0A]  sin[0ft])  sin[oA]) 


[149] 


[150] 


5.2  Miss  distance  termination  criterion 


The  estimated  missile-target  miss  distance  is  computed  based  on  the  assumption 
that  all  missile  and  target  accelerations  cease  from  the  time  of  estimation  to  the  point  of 
closest  approach.  The  relative  velocity  vector  is  then  projected  forward  to  define  the  point 
of  closest  approach  between  the  missile  and  the  target.  This  separation  is  the  projected 
miss  distance. 


From  geometrical  considerations,  the  miss  distance  dmi3ll  is  given  by 

dmiss  ~  RmT  m 

Vmt 


[151] 
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The  expression  for  Rmt%i  is  obtained  from  the  range  vector  in  the  LOS-fixed  coordinate 
system.  Since 

Rmt  =  Rmtxi  [152] 


then, 


•  — *  • 

Rmt%i  =  Vmt  —  Rmt&i  .  [153] 

Substituting  eq.  153  for  Rmt&i  and  eq.  85  for  Vmt  into  eq.  151,  the  resulting  expression 
for  the  miss  distance  is 


dmias  — 


[154] 


6.0  CONCLUDING  REMARKS 


DREV  is  currently  supporting  the  development  of  high-fidelity  simulation  mod¬ 
els  of  air-to-air  tactical  missiles  for  use  in  weapon  system  effectiveness  assessment.  For 
efficiency,  the  development  of  this  missile  simulation  capability  takes  advantage  of  already 
existing  model  components  when  available,  one  example  being  the  modeling  of  the  missile 
dynamics.  The  available  documented  model  of  the  missile  dynamics  uses  a  non-standard 
axes  system  where  x  is  forward,  z  is  lateral  and  y  is  upward.  This  model  was  modified  to 
adopt  the  standard  NED  axes  system  (x  forward,  y  lateral,  z  downward)  to  ensure  compat¬ 
ibility  with  existing  simulation  models  and,  in  particular,  with  the  DREV  library  of  missile 
components.  The  mathematical  relations  of  the  missile  dynamics  model  based  on  the  NED 
axes  system  have  been  developed  and  presented.  This  description  covered  the  Euler  trans¬ 
formations  between  the  six  coordinate  systems  required  for  the  missile-target  engagement 
simulation,  the  dynamics  of  the  rigid  airframe  and  of  the  seeker  (the  tracking  loop),  and 
the  computation  of  the  initial  conditions  and  of  the  miss  distance.  The  next  step  in  the 
development  of  the  required  simulation  capability  is  the  software  implementation  of  this 
missile  dynamics  model. 
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